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We analyse the recently available experimental data on direct photon productions from Au+Au 
collisions at ^snn=200 GeV RHIC and from Pb+Pb collisions at y/TNN—2.76 TeV LHC energies. 
The transverse momentum [pr) distributions have been evaluated with the assumption of an initial 
quark gluon plasma phase at temperatures ri=404 and 546 MeV with initial thermalisation times 
ri=0.2 and 0.1 fm/c respectively for RHIC and LHC energies. The theoretical evalutions agree 
reasonably well with the experimental observations. The thermal window for the LHC energy is 
justified through the ratio of the pr spectra of thermal photons to dileptons. 

PACS numbers: 25.75.-q,25.75.Dw,24.85.+p 



It is by now conventional wisdom, that collisions 
between two nuclei at RHIC as well as LHC en- 
ergies will lead to the formation of quark gluon 
plasma (QGP). Although the exact nature of QGP- 
hadron phase transition is still plagued by uncer- 
tainties, but there are several evidences that QGP 
is formed [J-Q at RHIC energies. It has long been 
recognised that promising signals of QGP are pho- 
tons and dileptons . The very nature of electro- 
magnetic interactions ensures that the thermomet- 
ric signals (7, ^jl^ , e+e") escape the QGP medium 
without any significant interaction, thus retaining 
the pristine information of QGP immediately after 
its formation. It has been pointed out by one of 
the authors and his collaborators, consistently for 
a long time [8l-[l^|. that the photons emanating 
from thermal medium show up in the pt window of 
1.5 < PT{GeV) < 3.5. It has also been pointed out 
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16j during the early days of QGP physics and 
now, that the ratio 7//i^/i~ eliminates some of the 
uncertainties associated with the input parameters 
used in the model. 

The concept of ratio has been revisited again in 
ref [13, [Si where the sensitivities of the of the ini- 
tial parameters have been studied extensively. The 
formation of a thermal medium is manifested by the 
flattering of the ratio beyond certain pT- Now that 
the direct photon data has come out of LHC we tend 
to analyse the data and test the above prediction at 
Viw7V=2.76 TeV. 

It is quite remarkable that the original theoreti- 
cal prediction p"5l - [l8j turns out to be valid also at 
LHC energies with added advantage that the very 
high energy of LHC turns the invariant mass M, 
of dileptons less relevant compare to (say) at RHIC 
energies. 

It should also be noted that there is already a 
tantalising hint of QGP formation even at SPS en- 
ergy il^ i2fli| . In the following work we present our 



findings of this interesting phenomenon. Thus over a 
wide range of energies going through orders of mag- 
nitude, it is quite exciting to note that pr window for 
thermal photons and dileptons remains the same as 
predicted earlier, only the range of window increases 
almost upto 4 GeV or more at LHC energy. 

In this present work we have evaluated the direct 
photon productions from (i) Au+Au collisions at 
^i^=200 GeV for 0-5% centrality and (ii) Pb+Pb 
colHsions at ^/snn=2.7Q TeV for 0-40% central- 
ity and finally compared with the recently available 
data [Mlfor RHIC and 24] for LHC energies. We 



analyse the data using our earlier approach, which 
has been used consistently to explain the SPS 
and RHIC data for different centralities. The 
invariant thermal yield of photons is evaluated from 
the production rate as in ref [13, UHl and using the 
(2+1) dimensional relativistic hydrodynamics. The 
photon contributions from prompt productions have 
been evaluated using the parton distribution func- 
tion as parametrised in Cteq6M. Different thermal 
sources of photon productions are considered as in 
ref 17. |25| and sources of prompt productions as in 
ref i3:l2al- The authors in 27] also mention the 
contribution of jet-plasma interaction to the high pT 
photons. Here the data seems to be explained nicely 
without considering the above phenomenon. The 
dilepton productions are evaluated as in [20| using 
the equation of state mentioned in sec-II. Then we 
evaluate ratio of the transverse momentum spectra 
of thermal photons to dileptons to seek out the ther- 
mal window at ^/snn=2.76 TeV. 

In section II the prompt and thermal photon pro- 
ductions have been dscussed. In the following sec- 
tion the expansion dynamics of the produced system, 
along with the initial conditions, is described. Fi- 
nally, the results have been discussed and compared 
with the available data. The summary is presented 
last. 
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On the basis of sources of production, photons 
emanating from the relativistic heavy ion colhsion, 
can be broadly categorized as follows; (i) prompt 
photons resulting from the interactions of the par- 
tons of the colliding nuclei, (ii) pre-equilibrium pho- 
tons, emitted before the medium gets thermalised, 

(iii) thermal photons originating from the interac- 
tion of thermal partons as well as thermal hadrons, 

(iv) photons produced the from the passage of jets 
through plasma and (v) the photons from the decay 
of long lived hadrons (tt^ — >■ 77, ij 77, etc). The 
aim of this work is to analyse the photon data from 
RHIC and LHC energies, where photons from the 
hadronic decays, have already been excluded. Out 
of the remaining categories, the pre-equilibrium con- 
tributions is negligible since thermalisation times(ri) 
for RHIC and LHC are very small. The photons 
due to the passage of jets through the medium is 
not considered here. In this work, we consider the 
prompt and thermal photons-the two main sources 
of direct photons, in terms of their transverse mo- 
mentum (pt) distribution. 

The prompt photons originate from the initial 
hard scatterings of the partons, primarily because 
of the compton processes {q{q) + — > q{q) + 7), 
annihilation process (5 + ^ — > 5 + 7) and quark frag- 
mentation {q{q) — > 9(9) + 7) following the scattering 
of partons of the nucleus in the colliding nuclei. In a 
complete and consistent NLOpQCD [0(q!Q;^)] calcu- 
lation, it is important to account for the prompt pho- 
tons arising from various 2 — )• 3 processes |28i, |2!|| . 
This calculation is comparable with the leading or- 
der calculation [l^ with a higher K factor. We use 
this calculation and scale it up by the number of 
binarry collisions for Au+Au and Pb-t-Pb collisions 
at y/sWN=200 GeV (RHIC) and 2.76 TeV(LHC) en- 
ergies respectively. The intrinsic fcy smearing is ig- 
nored here. 

The Pt distribution of thermal photons is ex- 
pressed in terms of invariant yield as follows; 
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N~f is the number of photons produced, y is the 
rapidity, where i represents QGP and Hadronic 

phases. ^ ipp^dy ) static rate of photon pro- 

duction at a temperature T from phase i. d^x is 
the four volume element and its evolution is taken 
care by relativistic ideal hydrodynamics. The pho- 
ton emission rate is given as follows; 

d^R gi"' 
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(see P-Q for details). ImH^^ is the imaginary part 
of the retarded photon self energy. Jb e is the ther- 
mal distribution function. (5^''=-Epoiarization ^''e'')- 
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FIG. 1: The pr spectra of direct photons from Au-|-Au 
collisions at yijvjv=200 GeV RHIC energy from differ- 
ent sources. 



The static emission rate of photons for QGP and 
hadronic phases are calculated from various partonic 
and hadronic processes. 

The photon emission rate from QGP due to comp- 
ton {q{q)g — > 9(9)7) and annihilation {qq — > 97) 
processes were evaluated [13, HH by using hard ther- 
mal loop (HTL) approximation [32|. Further it was 
found [33| that photon productions from the reac- 
tions, gq ^- gq-f, qq qqj, qqq qj and gqq gj 
contribute in the same order as annihilation and 
compton processes do. The suppression due to mul- 
tiple scattering durin g th e emission process which 
was not mentioned in |30ll3TI.[33t was later discussed 
in Ref. [sj]. The complete calculation of photon 
emission rate from QGP to order O(aas) has been 
completed by resuming ladder diagrams in the effec- 
tive theory [35|, which has been used in the present 
work. The emission rates for various processes are 
available in parametrised form in Ref. 136] . Here we 
use the temperature dependence of the strong cou- 
pling constant from Ref. [37[. 

An exhaustive set of hadronic interactions have 
been considered to evaluate the photon emission 
rate from a hadronic phase. The relevant processes 
are basically; (i) tttt — > p7, (ii) 7rp — > tt^ (with 
TT, p, uj, 4> and ai in the intermediate state [38[), 
(iii)7r7r rjj and (iv) nrj — nj. Also the radia- 
tive decay of higher resonance states [38^-'4l'| such 
as, p — >■ TTTT 7 and a; — > 7r7 produces photons. The 
corresponding vertices can be obtained from vari- 
ous phenomenological Lagrangians as described in 
detail in Ref. (38.-41i] . The contributions from the 
reactions involving strange mesons like ttK* — > Kj, 
ttK — )• K*j, pK — >• and KK* — )• Kj, have been 
pointed out in [i^l Contributions from other decays. 
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such as iC* (892) K-i,(j) 777, 6i(1235) 7r7, 
02(1320) ^- 717 and ^1(1270) 7r7 have been 
found to be small [i^ for pT > I GeV. We con- 
sider all the processes for our calculation. All the 
isospin combinations for the above reactions and de- 
cays have properly been taken into account. The 
effects of hadronic form factors(dipole) [l^] have 
been incorporated. 

When two energetic heavy nuclei collide with each 
other, a large amount of energy and thus entropy is 
channeled into a small volume. The produced sys- 
tem achieve thermalisation in a short time because of 
the secondary interactions and the matter then ex- 
pands due to the high internal pressure. The space- 
time evolution of the matter has been studied using 
ideal relativistic hydrodynamics ^44*] with longitudi- 
nal boost invariance (45|] and cylindrical symmetry. 
In the present analysis of Pb-|-Pb collisions we con- 
sider the initial energy density (e(Ti,r)) and radial 
velocity i^i ji^yy) p rofiles similar to our earlier stud- 
ies [ij, [13, [20, lla) . Transition temperature which is 
an input parameter to the hydrodynamic calculation 
is taken as Tc=163 MeV according the recent works 
of lattice QCD computation by the HotQCD Collab- 
oration [46|] and by the Bielefield Group [4^. We 
also consider Tc= 170 MeV as predicted by the lat- 
tice computation from other group '^s'] to compare 
the results. The initial temperature Ti and initial 
thermalisation time Ti are constrained from the fol- 
lowing equation poj : 
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where C(3) denotes the Riemann zeta function, R 
is the radius l.l{Npart/'2,y^^ : Npart is the num- 
ber of participant nucleons ] of the colliding sys- 
tem , Ti is the initial thermalisation time and 5e// 
is the statistical degeneracy taken to be 37, con- 
sidering a 2-flavor QGP to be produced initially. 
dN / dy{b) is the hadron (mostly pion) multiplicity 
for a given impact parameter h, which is obtained 
from the Glauber model for 0-40% collision central- 
ity of Pb-I-Pb colHsions at ys]^=2.76 TeV. The 
initial temperature obtained from the above equa- 
tion Eql3]for dN/dy{b)=1275 and t,;=0.1 fm is 546 
MeV. 

The heavy ion collision experiment at LHC en- 
ergy maps the region , where the value of baryonic 
chemical potential, /is is small and the temperature 
is rather high. Guided by the lattice computation, 
which predicts a crossover from hadronic to par- 
tonic phase for ^^=0, we construct an equation of 
state (EOS) for the evolution dynamics which gives 
crossover like transition. We consider the entropy 
(sq) of the QGP phase (for T > T^) using BAG 
model EOS and the entropy {sh) of the hadron gas 
(considering non-interacting hadrons and their res- 
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FIG. 2: Direct photons measured by PHENIX collabo- 
ration for 0-20% centrality at ^5]^=200 GeV 



onances up to mass ~ 2.5 GeV for T < Tc). The 
entropy during the transition region is parametrised 
using a tan-hyperbolic function [50| as follows. 



and fqiT) = i(l + tanh^-^) 



(4) 



Here T is the width parameter and assumes a finite 
value for crosssover transition. This value can be 
tuned to zero to the first order transition. Here the 
width parameter is taken to be r=20 MeV. EOS 
which is an input to the hydrodynamic calculation 
has been constructed for a transition temperature 
rc=163 MeV and width parameter r=20 MeV. The 
effective degeneracy g(T) of the system which is ob- 
tained from the entropy s{T) (Eq. |4]). 

The Pt spectra of direct photons have been evalu- 
ated for Au-f Au collisions at y/sWN^2QQ GeV RHIC 
energy. For 0-5% of the centrality of the collision, 
the total pion multiplicity is taken as dN / dy=llQQ. 
To evaluate the thermal contribution, we consider 
an initial temperature Ti=404 MeV, which is ob- 
tained from the Eq. [3] for an initial thermalisation 
time Ti=0.2 fm. An equal time freeze out scenario 
for all the hadrons is assumed and the value of 
freeze out temperature Tf is taken to be 120 MeV, 
which is contrasted from the pion and kaon spec- 
tra [l3| . The prompt photon contributions have 
been evaluated from the annihilation, compton and 
quark-fragmaneation processes using pQCD calcula- 
tions [l3|- The calculation is done for A^co/;=910 (ob- 
tain from Glauber model calculation for 0-5% cen- 
trality) and tTi„=42 mb. We use the parametrised 
parton distribution function from CteqdM . To ac- 
count for for the next-to-next leading order contri- 



bution, the value of and Ki 



brem 



is taken to be 1.4 
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FIG. 3: The pion spectra from Pb+Pb collisions at 
y/sNN=2.76 TeV LHC energy fitted with blast wave with 
temperature T/^f f3 MeV and ^=0.9 



and 2.2 respectively for RHIC energy. The results 
are shown in Fig. [T] The solid cirlces represent the 
recently published direct photon data from Au+Au 
collisions (0-5% centrality) by PHENIX collabora- 
tion [21]. The long-dashed and dotted lines repre- 
sents the pQCD and thermal contributions respec- 
tively. The solid is for the sum total of both. T^, is 
taken to be 163 MeV. 

The Fig. [2] displays the direct photon productions 
for 0-20% centrality of Au-f Au colhsions at RHIC 
energy [22| as per our earlier evaluation [13|]. The 
dotted line represents the thermal contribution for 
lattice EOS ^ with T,=300 MeV and t,=0.5 fm. 
The prompt photons are evaluated by scaling the 
NLO calculation of Gordon and Vogelssang's [26| 
prediction for p-p collisions with the number of colli- 
sions for Au-|-Au. The solid line represents the total 
photons. 

While evaluating the px distributions for direct 
photons at LHC energy, same time freezeout sce- 
nario for all the hadrons is assumed with Tf— 113 
MeV, which is constrained by the pion data. Fig. [3] 
displays the blast wave explanation of the pion data 
from Pb+Pb coUisions at ^/snn— 2.76 TeV (For data 
see (111 ) with T=113 MeV. 

In Fig. 2] the direct photon data from Pb-|-Pb colli- 
sions for 0-40% centrality at y/sNN=2.76 TeV mea- 
sured by ALICE collaboration (solid circles) [l^l 
is shown. The dotted line represents the pQCD 
calculation of annihilation and compton processes 
(prompt) using CteqSm for A^coi;=853(for 0-40%). 
The solid line with open circles show the contribu- 
tions from quark-fragmentation to the prompt pho- 
tons. The long-dashed line is the sum over all the 
hard processes. Kj=3.2 for annihilation and comp- 
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FIG. 4: The pt spectra of direct photons from Pb-|-Pb 
collisions at ^snn='2-76 TeV LHG energy from different 
sources. 
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FIG. 5: The pr spectra of direct photons from Pb+Pb 
collisions at ^snn=2.76 TeV LHC energy. The solid 
cirlces are the data points measured by ALICE collab- 
orations f2ll|. The solid line represents the theoretical 
evaluation. 



ton processes and -ft'f,rem=2.8 for fragmentation pro- 
cess are considered to take the NNLO contributions 
into account. Considering Ti to be 546 MeV, 
which is obtained from Eq. |3] for 0-40% centrality 
and T,=0.1 fm, T/=113 MeV and Tc=163 MeV, the 
thermal contributions are evaluated from QGP and 
hadrongas. This is displayed in Fig. |3](solidline). As 
depicted from Fig. SI the thermal contribution domi- 
nates the invariant yield upto p-p ~ 4 GeV. The spec- 
trum beyond pr > 4 GeV is dominated by pQCD 
calculations clearly distinguishing thermal window 
1.5 < priGeV) < 4.0. The solid line in Fig. [5]dis- 
plays the sum of thermal and pQCD contributions 
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explaining the data within the ambit of considered 
parameters. 

It is important to mention the sensitivity of the 
results due to the initial parameters. The transition 
temperature Tc=163 MeV is restricted from the re- 
sults of lattice QCD computation by HotQCD and 
Bielefield group G^, H^]. However the computa- 
tion as given in [48'| predicts the value around ~ 
170 MeV. To differentiate the effect of we evalu- 
ate the photon spectra for LHC energy with Tc— 170 
MeV and found very little change in the pT spec- 
tra. At pT=2 GeV, the resuhs for T^^UO MeV is 
- 1.2% higher than rc=163 MeV. The sensitivity 
of Tc is also found to be small in our earlier work 
[isj . The sensitivity of results to the initial tem- 
perature (in GeV) and thermalisation time (in fm) 
(Ti, Ti) together is shown in Fig. [5] Comparing the 
multiplicity we consider (T^ = 433, = 0.2) and 
reevaluate the photon spectra. The solid line is for 
{Ti — 546, Ti = 0.1) and the long-dashed line is for 
(T, = 0.433, T, = 0.2). At pt^2 GcV the result 
differs by 7.54%. But the results are found to be 
sensitve to Tf. In Fig. [7] the prediction for direct 
photons at ■y/SArAr=5.5 TeV is given with an ini- 
tial QGP temperature ri=843 MeV and ri=0.08 fm. 
The charged particle multiplicity is taken to be 2500 
for 0-5% centrality. It shows the dominance of ther- 
mal photons up to pt=4 GeV of the direct photon 
spectra. 

It has been argued since long [H, [l^ , that the 
ratio of the pr spectra of thermal photons to lep- 
ton pairs forms a plateau in the pt window 1.5 < 
pxiGeV) < 3.5 for different mass bins of lepton pairs 
and it is less sensitive to the input parameters. In 
ref [l3| the ratio has been revisited and the sensitiv- 
ities to different input parameters have been studied 
extensively. Here we study the ratio for ^/snn =2.76 
TeV LHC energy for different mass window of the 
lepton pairs and plot the results in Fig. [Sj The ther- 
mal dilepton yield is calculated using the approach 
mentioned in [2Cl [s^ , that explains the dimuon 
data from In-fin collisions measured by NA60 col- 
laboration. A plateau is observed beyond pr ^ 1.5 
GeV in Fig. 13 

The plateau in the ratio appears because of the 
following reason. The emission rate of photon and 
dilepton production show their momentum depen- 
dencies through their thermal phase space factor 
fBE{E,T) (^ exp{- E/T)). Where E = Mrcoshy, 
Mt = y/pT + M'^ & y = tanh-^p^/E. p^ and 
y take their usual meanings. The nature of the 
spectra depends on M . When we consider the ra- 
tio of their momentum spectra at high p^ (px » 
Mj&zMt = Pt), for a static system, then the real 
photon {M^ — 0) and dilepton (M^ 7^ 0) have sim- 
ilar momentum dependence. Thus a plateau is ob- 
served in ratio for the static system. However for 
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FIG. 6: Sensitivity in the production of thermal photons 
to Ti for LHC energy. The long dashed line represents 
the results for Ti=433 MeV and fm. Solid for 

Ti=546 MeV and Ti=0.1 fm 
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FIG. 7: Direct photon procutions at ^sjvjv=5.5 TeV 
LHC energy. 



a realistic case of an expanding system the ther- 
mal phase space factor goes as /b_e ~ exp{—u^p'^), 
where and are the 4- velocity and -momentum. 
The spectra, then, depend on the radial flow Vr along 
with M. It has been observed in [l^ that for (i) 
Pt » M, the radial flow of photons and dilep- 
tons are similar, and thus the plateau is achieved, 
(ii) If the large M pairs originate from early time, 
when the flow is small, the ratio which includes the 
space time dynamics will also be close to a static 
case showing a plateau, (iii) But the plateau dissa- 
pears when (pT ~ M) or radial flow is large. Here 
we present the ratio for different lepton pair mass 
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FIG. 8: Ratio of the invariant yield of thermal photons to 
dileptons at -^sjv]v=2.76 TeV LHC energy for different 
invariant mass windows. 



windows for the range of 0.3 < M{GeV) < 1.3 at 
Vs]^=2.76 TeV. 

It may be noted that the ratio of photons to dilep- 
tons obtained from pQCD calculation (non-thermal 
origin) doest not show a plateau, which is analysed 
extensively in 17 1. 



In summary, we say, that the measured direct pho- 
ton spectra from Au-|-Au coUisions at ^/snn— 200 
GeV for 0-5% centrality and from Pb-|-Pb colli- 
sions at y/sNN=2.76 TeV for 0-40% centrality have 
been reproduced by ideal relativistic hydrodynamics 
with an EOS which represents crossover like transi- 
tion. Initial QGP phase is assumed for both RHIC 
(T,=404 MeV, n=0.2 fm) and LHC (r,=546 MeV, 
r,=0.1 fm). 

It is observed somewhat remarkably that over a 
wide range of energies going through orders of mag- 
nitude, from (even SPS) RHIC to LHC, the ther- 
mal photons populate the same pt window 1.5 < 
pxiGeV) < 3. The ratio j/fi'^ii~ for the invariant 
mass widow 0.3 < M{GeV) < 1.3 and pt window 
1.5 < pxiGeV) < 3, turns to a plateau indicating 
the onset of QGP in the heavy ion collision. The 
formation of thermal window is manifested in the 
ratio within 1.0 < priGeV) < 4.0 at ysKw=2.76 
TeV LHC energy. 
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